function output = my_getThirdSeries(MATRIX,LAT,LON)
% output = my_getthirdseries(MATRIX,LAT,LON)
%
% --------自定义函数--------
% 函数名称：my_getThirdSeries
% 函数说明：获取lat-lon(-time)矩阵中指定点（lat，lon）的时间序列
% 函数输入：MATRIX,LAT,LON
%   MATRIX： lat-lon-time矩阵
%   LAT：   纬度，或者行号，可以是列向量。
%   LON：   经度，或者列号，可以是列向量。
% 函数输出：output
%   output： 列向量或列向量（时间序列）构成的矩阵
% 作者：sid
% LOG：
%    2021/10/26 ver 1.0: create
%    2021/10/26 ver 2.0: update 支持在lat-lon矩阵上取样。
switch ndims(MATRIX)
    case 2
        LinearIndex = sub2ind(size(MATRIX),LAT,LON);
        output = MATRIX(LinearIndex);
    case 3
        Matrix = my_shiftdim(MATRIX,1);
        LinearIndex = sub2ind(size(MATRIX,1,2),LAT,LON);
        output = Matrix(:,LinearIndex);
end
end